Flat bands, Dirac cones and atom dynamics in an optical lattice 
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We study atoms trapped with a harmonic confinement in an optical lattice characterized by a 
flat band and Dirac cones. We show that such an optical lattice can be constructed which can be 
accurately described with the tight binding or Hubbard models. In the case of fermions the release 
of the harmonic confinement removes fast atoms occupying the Dirac cones while those occupying 
the flat band remain immobile. Using exact diagonalization and dynamics we demonstrate that a 
similar strong occupation of the flat band does not happen in bosonic case and furthermore that 
the mean field model is not capable for describing the dynamics of the boson cloud. 
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I. INTRODUCTION 

Recent development of manipulating cold atoms in op- 
tical lattices has opened new avenues for studying corre- 
lation and band structure effects in a controlled way^E]- 
One, two and three-dimensional lattices with a variety 
of lattice structures have been proposed, including the 
kagome lattice[6-9 , which was originally suggested by 
Syozi[rO and studied intensively in the condensed mat- 
ter physics mainly due to its magnetic properties [1 IHl4] . 
In a simple tight binding (TB) model the kagome lat- 
tice is characterised by a completely flat band, meaning 
that the energy is independent of the wave vector k and 
that electrons occupying this band are localized. Flat 
band Hubbard model is a paradigm for ferromagnetism, 
recently reviewed by Tasaki [15 . 

Kagome lattice is one example of a whole class of 
two (2D) and three-dimensional (3D) lattice structures 
which, in the TB model, produce a flat band [TH [17]. In 
2D these structures include square and hexagonal lat- 
tices and can have one or several flat bands as well as 
crossing points where two bands open as circular cones 
leading to Dirac fermions, as in graphene[18 . The Dirac 
fermions are also known to lead to interesting magnetic 
properties [19]. Allowing also the p-states in each lat- 
tice site to be occupied even more freedom to create flat 
bands and Dirac cones is obtained. Wu et a/[8 have 
studied the honeycomb (graphene) optical lattice with 
p-electrons and demonstrated the atom localization due 
the flat band of the system. 

One of the simplest flat band lattices is the edge- 
centered square lattice with three atoms in a unit cell [20 , 
depicted in Fig. [ij In the TB model this lattice has three 
bands, the center band is flat and the upper and lower 
bands meet the flat band at the corners of the Brilloun 
zone where Dirac cones open. This means that at the 
same energy the particle can have an infinite effective 
mass (flat band) or zero effective mass (Dirac fermions). 
Recently Shen et al. showed that these massless Dirac 
fermions can exhibit perfect all-angle Klein tunneling 
[2T] . This is the lattice of interest in the present work. 

We will first show that laser field can be used to create 
a lattice which accurately produces a flat band and Dirac 



cones meeting at the same energy. The band structure 
is then fold to a TB model and for interacting atoms to 
a Hubbard model. In the case of fermions with only 
one spin state (spinless fermions) the Hubbard model 
equals to the noninter acting fermions since the Pauli ex- 
clusion principle prevents two similar fermions to occupy 
the same lattice point and on-site interaction has no ef- 
fect. This provides us a simple way to study a dynamical 
problems, for example, what happens when a harmonic 
confinement, keeping the atoms in the central region of 
the lattice, is removed. 

In the case of bosons the many-particle problem of the 
Hubbard model becomes more complicated. In this case 
we solve the problem with the Hubbard model for a small 
system of only three particles. However, this already 
demonstrates that (i) the boson system behaves differ- 
ently than the fermion system and (ii) that frequently 
used mean field models (e.g. Gross-Pitaevskii), where the 
boson system is described with only one single particle 
wave function, fails to describe the dynamics correctly. 

The results show that, when the harmonic confinement 
is removed, part of the atoms fly away fast while part 
remain stuck in the immobile flat band states. We an- 
ticipate that this kind of experiment can be made in the 
near future using both bosonic and fermionic atoms. 

The edge-centered square lattice with three atoms in 
a unit cell is illustrated in Fig. [l] We consider fist a 
simple TB model with only one state per lattice site, 
assume only nearest neighbour hopping and neglect the 
so-called differential overlap between neighbours. The 
band structure can be easily solved [22] by diagonalizing 
a 3 X 3 matrix for each 2D wave vector q = (qx^Qy), 
resulting to energy levels 

ei(q) = —t^y4 + 2 cos{qxa) + 2 cos{qya) 

^2(q) = (1) 

e3(q) = ty^A + 2cos((7a,a) + 2cos((7^a), 

where t is the strength of the hopping integral between 
nearest neighbours and a the lattice constant. The corre- 
sponding wave vectors Vi = {vc, '^ei, '^e2), where the three 
components refer to the corner and edge sites in the unit 
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cell, are {s = V4 + e^^-^ + e"^*?-^ + e*^^/^ + e"*^^'^) 

Vi(q) = (-5, 1 + l + e-^*?-^) 

^2(q) = (0, l + e*^-'^, l + e^^^«) (2) 

The energy bands consist of a flat band at zero energy 
and two bands symmetrically below and above the flat 
band. All bands meet at points at the corners of the Bril- 
louin zone, = Qy = ir/a. The band structure is shown 
m Fig. [1] The wave function of the flat band has zero 
amplitude at the corner point of the square lattice. This 
is the reason of the flatness: The particles in these states 
do not hop to the neighboring sites and can be viewed 
as localized. The effective band mass of the particles is 
thus infinite. 
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FIG. 1: Up left: Edge-centered square lattice and its unit cell 
shown with a dashed line. Up right: Pptimized potential, Eq. 
(3), with equally deep minima shown with red color. Down 
left: Three lowest bands corresponding to the potential of 
Eq. (3). The cusps are at the corners of the Brillouin zone. 
Down right: Arrows show the suggested directions of the laser 
fields causing potentials like cos^x, cos^{x-\-y) and cos^(2x), 
corresponding to minima at the sites of the red, blue and 
green dots, respectively. 

Close to the zero energy the two other bands form 
"Dirac cones", and the velocity of the particles is in- 
dependent of the energy, ei^s{kx — Tv/a^ky — 7r/a) ~ 

^ta^k^ -\- ky/4: = where k is the distance from 

the point of the cusp. The linear dispersion relation 
means a zero effective mass. Interestingly, in this lat- 
tice we have at the same energy region particles with 
zero effective mass and infinite effective mass. Note that 
the cones appear to be symmetric to rather high energies 
indicating that most of the particles in the lowest and 
uppermost bands have nearly equal velocities. 

Next we will consider how such a simple TB lattice 
could be constructed for cold atoms in an optical lat- 
tice. We do not go in the details of the interactions 
between atoms and laser fields ^j-[3], but simply assume 



that the effective potential caused by the standing elec- 
tromagnetic waves can be described with trigonometric 
functions (cos^(x) or equivalently cos(2x)). Clearly a 
square lattice can be made as — cos^(ax/7r) — cos^(ay/7r), 
having minima at lattice sites R = (na^ma). Minima 
at the edge centers can be constructed with functions 
cos^(a(x + y)/2) + cos^(a(a:: — y)/2. However, with a lin- 
ear combination of these two it is not possible to con- 
struct proper potential barriers between corner end edge 
sites and an additional shorter wave length function of 
— cos^(2ax/7r) — cos'^ {2ay / ir) has to be added. The po- 
tential optimized to give the band structure similar to 
that of the TB band structure is 

Voi^{x,y) = ■\/o[-0.9(cos^[ax/7r] + cos^[a^/7r]) 

+0.496(cos2[a(x + y)/2] + cos^[a{x - y)/2]) 
+(cos^[2ax/7r] + cos^[2ay/7r])], 

(3) 

where Vq gives the energy scale (see below). The resulting 
potential is shown in Fig. [l] The lowest three bands of 
the band structure (computed using a real space mesh of 
26 X 26 points) are shown in Fig. [l] and they are nearly 
identical to those of the simple TB model. The parameter 
Vq is related to the mass of the atoms m and the lattice 
constant a as Vq ~ 497fi^/(2ma^) (where the numerical 
factor 497 comes from fitting of the band structure to 
that of the TB model). 

The overall width of the three lowest bands is 
7 .Qbh^ / (2mo?) while the width of the flat band is only 
0.11fi^/(2ma^), meaning that the dispersion of the flat 
band is only 1.5 % of the total width of the three bands. 
The potential of Eq. (3) has, of course, an infinite number 
of other bands, but there is a huge gap of 140fi^/(2ma^) 
between the three lowest bands and the next band. With 
this example we have demonstrated that a proper combi- 
nation of cosine potentials can produce a band structure 
which, with a high accuracy, gives similar energy bands 
as the simple TB model. Figure [T] shows schematically 
the laser arrangement needed for making such an optical 
lattice. 

We assume that the atoms trapped by the optical lat- 
tice contact with a repulsive short range potential and 
assume the interaction is so strong that only one atom 
is allowed in each lattice site. The many-particle prob- 
lem then reduces to the Hubbard Hamiltonian. We first 
consider fermionic atoms with only one spin state (spin- 
less fermions). In this case the many-particle problem 
is noninteracting, since the Pauli exclusion principle al- 
ready prevents two atoms to occupy the same site. 

In the actual experiments the optical lattice is super- 
imposed on a harmonic confinement which localizes the 
particles (with repulsive interaction) in the central re- 
gion of the lattice. Putting enough atoms in the lattice 
the harmonic confinement guarantees that in the central 
region the occupation of each lattice site is nearly one. 
This means that all the three bands (which are bend in 
the harmonic confinement) are filled at the central region 
of the harmonic confinement. 

Our interest is to see what happens when the harmonic 
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confinement is removed but tiie optical lattice is kept in 
place. This dynamical problem is easily solved for N 
spinless fermions. The procedure is as follows: We solve 
the TB model in a large finite lattice with a harmonic 
confinement and fill the lowest N single particle states. 
This is the fermionic ground state ^o(O) at time zero. 
Next we remove the harmonic confinement and expand 
the states with the confinement in terms of the states 
calculated without the confinement. The time depen- 
dence follows then from the time dependences of each 
individual state as 

Mt) = T.^je-'^''^''^°{0), (4) 

j 

where Ej is the energy of the many-particle state j. In 
the case of spinless fermions where each state is a Slater 
determinant of single particle states the time dependence 
further reduces to the time dependences of the single par- 
ticle states. 
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(c) t = 5 (d) t - 7.5 
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FIG. 2: Time evolution of 50 fermionic atoms when the ar- 
monic confinement is removed, (a) shows the initial state 
with the harmoic confinement of strength Uh — OA (in units 
of fi? /2ma^)^ and (d) the 'final' state when atoms trapped 
by the flat band stay localized while other otoms have flown 
away. 

In practice we solved up to 100 atoms in a lattice 
of 25 X 25 unit cells. We use TB units with the hop- 
ping parameter t = 1 and the nearest neighbour distance 



6 = a/2 = 1. The harmonic confinement adds to diag- 
onal terms to the TB Hamiltonian matrix: Placing the 
center of the finite lattice at the origin, the harmonic 
confinement changes the energy level of a lattice site by 
/7hr^, where is the strength of the harmonic confine- 
ment and r the distance of the lattice site from the origin 
(in units of the lattice constant). Figure [2] shows an ex- 
ample of such a simulation. The Harmonic confinement 
localizes the atom cloud at the center of the finite lattice. 
The maximum density is at the edge sites of the unit cell 
with occupation about 0.8, while at the corner sites it is 
about 0.6. When the harmonic confinement is removed, 
the cloud fast expands. Fig. |2] (b) and (c), but part of 
the atoms stay trapped at the edge sites, (d), forming a 
lattice of localized atoms. The occupation of the edge 
states remains close to 0.5. 

The explanation of the observed dynamics is a direct 
result of the single-particle properties of noninteracting 
fermions: Those occupying the bands 1 and 3 are mobile 
and fly fast away (at the dirac cones even with the same 
velocity), while atoms occupying the flat band are im- 
mobile and stay in place. Notice that the corner states 
of the unit cell are emptied since the flat band does not 
have any amplitude in these sites. 

The bosonic many-particle problem is usually simpler 
than the fermionic case due to the symmetry of the wave 
function. In our case, however, this is not the case. Even 
if we assume an infinitely strong contact interaction the 
bosonic case remains a true many-particle problem, al- 
though the boson and fermion problems are related |23]. 
A common approach is to use a mean field model for 
the Hubbard Hamiltonian (this is closely related to the 
Gross-Pitaevskii [24| [25] model): The interaction term 
U{ni — l)ni is replaced with U{{ni) — l)n^ (where fii 
is the occupation number operator). The first approxi- 
mation is to assume that all the bosons occupy the same 
quantum state determined by the external confinement. 

Taking the interaction U very large, we prevent the 
occupation of any lattice site to be much larger than 
one. In the Harmonic confinement we can then get ini- 
tial atom density to be qualitatively similar to that of 
the fermionic case and one would expect that expanding 
the self-consistent bosonic ground state in terms of the 
single particle band states also the flat band is markedly 
occupied. However, this is not the case. The occupa- 
tion of the flat band state remains less than 2 % even 
when the interaction strength U is extrapolated to infin- 
ity. This means that when the confinement was removed, 
practically all the atoms of the cloud escape fast. 

In order to understand if the boson system behaves 
truly differently than the fermion system, or if it is the 
mean field approach which is not appropriate for our sys- 
tem, we performed a small scale exact diagonalization of 
the Hubbard Hamiltonian. The test size was 3x3 unit 
cells with periodic boundary condition, i.e. 27 sites alto- 
gether as shown in Fig. |3] We assumed infinitely strong 
repulsive interaction between the atoms {U = oo) which 
reduced the boson basis to states having at most one 
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FIG. 3: Left: Finite lattice used, wiht periodic boundary con- 
ditions, for solving the Hubbard model for three bosons or 
fermions. The initial potential at the cites denoted by red and 
blue were — 20t and — 5t, respectively. Right: Time dependece 
of the occupation of red and blue sites after the potential of 
these sites were put to zero. Solid lines are results for fermions 
and the thick dashed lines for bosons. The thin dashed lines 
are results for a single particle with the same initial density as 
the exact boson result, indicating that the mean field solution 
fails in this case. The time is in units oi h/t 

atom in a lattice site. Intially the atoms were localized 
in a small region with strong attractive potentials in five 
adjacent sites as shown in Fig. |3] By diagonalizing the 
Hubbard Hamiltonian for three bosons (and fermions) 
with and without the attractive potentials at five sites, 
the time evolution could be solved exactly using Eq. (4) . 

The time evolution is shown in Fig. [3] for bosons and 
fermions. In the case of fermions the results is qualita- 
tively similar to that of larger systems described above. 
After a short time the atom density at the corner site 
(red) decreases to nearly zero while at the edge sites 
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